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We have numerically studied the thermodynamic properties of the spin | XXZ chain in the pres- 
ence of a transverse (non commuting) magnetic field. The thermal, field dependence of specific heat 
and correlation functions for chains up to 20 sites have been calculated. The area where the specific 
heat decays exponentially is considered as a measure of the energy gap. We have also obtained the 
exchange interaction between chains in a bulk material using the random phase approximation and 
derived the phase diagram of the three dimensional material with this approximation. The behav- 
ior of the structure factor at different momenta verifies the antiferromagnetic long range order in 
^/-direction for the three dimensional case. Moreover, we have concluded that the Low Temperature 
Lanczos results [M. Aichhorn et al., Phys. Rev. B 67, 161103(R) (2003)] are more accurate for 
low temperatures and closer to the full diagonalization ones than the results of Finite Temperature 
Lanczos Method [J. Jaklic and P. Prelovsek, Phys. Rev. B 49, 5065 (1994)]. 

PACS numbers: 75.10.Jm, 75.40.-s, 75.40.Cx, 75.40.Mg 



I. INTRODUCTION 

Quantum phase transition in strongly correlated elec- 
tron systems which is the qualitative change in the 
ground state properties versus the a parameter in the 
Hamiltonian (a magnetic field, amount of disorder, • • • ) 
have been at the focus of research recently^'^. Spe- 
cially, field induced effects in the low dimensional quan- 
tum spin models have been attracting much interest 
from theoretical and experimental point of view in re- 
cent years^i^i^i^iii^i^ii^iiiii^. Thc magnetic properties of 
a system with axial anisotropy depends on the direction 
of the applied field. For instance, the magnetic properties 
of the antiferromagnetic one dimensional spin 1/2 XXZ 
chain in the longitudinal field is quite different from the 
case of transverse field. The longitudinal field commutes 
with the rest of Hamiltonian and preserves the integrabil- 
ity of the model by Bethe ansatz while a transverse field 
does not commute and the model is no longer integrable. 
The transverse field induces the antiferromagnetic long 
range order in perpendicular direction to the field and 
causes a quantum phase transition to the paramagnetic 
phase at the critical point. Experimental observations^^ 
justify this effect in the quasi-one dimensional compound, 
CS2C0CI4. 

The qualitative change in the ground state at the quan- 
tum critical point (QCP) for zero temperature (T = 0) 
affects the finite temperature properties of the model 
close to the QCP. The quantum critical properties of the 
Ising model in transverse field has been extensively stud- 
ied with different approaches^' However, there is no 
work on the finite temperature properties of the XXZ 
model in the presence of a transverse field. It is our aim 
to study the thermodynamic behavior of the anisotropic 
Heisenberg chain in a transverse field which is the model 



Hamiltonian for the mentioned behavior. The Hamilto- 
nian for this system can be written as 

N 

H = jY^is^s-^, + s\s\^^ + Asfsf+i - hst), (1) 

1=1 

where J > is the exchange coupling, < A < 1 is the 
anisotropy in z-direction, h is proportional to the trans- 
verse magnetic field and is the a-component of Pauli 
matrices at site i. Experiments^^ showed that CS2C0CI4 
is an realization of this model with J = 0.23 meV and 
A = 0.25. 

We will shortly discuss the different approaches to fi- 
nite temperature properties of the lattice model using 
Lanczos method in the next section. We then implement 
an appropriate approach to get the thermodynamic prop- 
erties of the XXZ model in the transverse field. In this 
paper, we have studied the finite temperature properties 
of a chain with = 20 sites in magnetic fields h = . . .5 
and with anisotropics A = and A = 0.25. 



II. FINITE TEMPERATURE LANCZOS 
METHOD 

The Lanczos diagonalization method is a powerful nu- 
merical tool to study the properties of the quantum many 
body systems on finite clusters. It is usually used to get 
the ground state of a quantum model with high accu- 
racy. However, this technique can be extended to finite 
temper at urei^ii^, where one can study the energy spec- 
trum of the system and its behavior under changing sys- 
tem parameters such as anisotropy and magnetic field. 

In the Lanczos algorithm, the ground state can be ob- 
tained with very high accuracy. However, to get finite 



temperature behavior we need to take average over the 
whole Hilbert space (L) with Boltzmann weights, i.e. 
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It has been proposed that the Lanczos procedure can 
be used to get rehable approximation for the finite tem- 
perature (T 7^ 0) properties of the lattice models^^. In 
this approach the Hilbert space is partially spanned by 
different random initial vectors for several Lanczos proce- 
dures. This approach is called Finite Temperature Lanc- 
zos Method (FTLM)^^ which is based on the following 
equations. 

R M 

=^ ^EE^"'"'" ^''i*™ )(*™ i^i'-) (3) 
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where P = l/ksT (ks = 1) is the Boltzmann constant 
and ~ 10 is the total number of random samples 
with different initial vectors |r). The eigenvectors and 
the corresponding eigenvalues of the tridiagonal matrix 
are \^^) and respectively for m = 1, . . . , M. The 
FTLM procedure works well for finite temperatures, how- 
ever, it does not converge to the ground state expectation 
value ((V^olOIV^o)) as T ^ because of statistical fiuctu- 
ations. i.e. 



{O) = X^(*o|0|r)(r|vl/ o)/E(*o|r)(r|vl/ o) (5) 



This problem can be solved by a symmetric algorithm 
which is explained below. 



A. Low Temperature Lanczos Method 

The Low Temperature Lanczos Method (LTLM) was 
proposed as a symmetric algorithm to remedy failure of 
FTLM at low temperatures^^. However, LTLM needs 
more CPU time and memory but requires fewer steps in 
low temperatures^. LTLM can be formulated similar to 
the case of FLTM with the following equation. 
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(6) 



where the partition function {Z) is calculated in the same 
way as in Eq. (|3|). Using this method, as T ^ we have. 
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FIG. 1: The x-component spin structure factor at momentum 
TT versus temperature for a chain of length = 10, 12 using 
FTLM, LTLM and exact diagonalization (ED). (Color online) 



(O) = £(*o|r>(r|*o)(*o|0|*o)/E(*ol^)<^l*o) (7) 
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which is equal to the ground state expectation value. 
It is thus more accurate to get the low temperature quan- 
tum properties. 

III. FINITE TEMPERATURE RESULTS 

A. Comparison between LTLM and FTLM 

We are interested in the finite temperature properties 
of the Hamiltonian defined in Eq. ([1]) and specially in the 
quantum critical behavior as T ^ 0. We have thus imple- 
mented the LTLM method to get the finite temperature 
properties which is accurate close to T = 0, although it 
consumes more CPU time and memory. To justify this, 
we have plotted the x-component spin structure factor at 
momentum tt (G^^(7r), see Eq. ([9|)) versus T in Fig. [TJ In 
this figure, we have shown both the FTLM and LTLM 
results for a chain of length N = 10, 12 together with 
the exact diagonalization (ED) results for comparison. 
The discrepancy of FTLM in low temperature (T ^ 0.2) 
in comparison with the ED results is clear specially for 
TV = 12 as discussed in previous section. Fig. [1] shows 
that the LTLM results converge to the ED ones as T ^ 
and are in good agreement for higher temperatures. 

B. Specific Heat 

We have plotted the specific heat of the XXZ chain in 
the transverse field (h) versus temperature (T) in Fig. [2l 
The specific heat is computed using 



Cy 



J 



HE') - (Ef). 



(8) 



3 



The data comes from the Lanczos algorithm with ran- 
dom samphng using R = 100, M = 30 — 100 (to get 8 
digits accuracy in the first excited state energy) and for a 
chain of length TV = 20 with A = 0.25. For h < he ^ 3.3 
the specific heat shows exponential decay as T goes to 
zero. This is in agreement with the presence of a finite 
energy gap in this region. The model is gapless at = 0, 
however, our data on a finite size does not show this be- 
havior because of the finite size effects. The finite size 
effects also appear as some level crossing between the 
first excited state and the ground state in finite system 
which become degenerate ground state in the thermody- 
namic limit {N oo). This level crossing causes small 
oscillations close to T = and for h < he. The finite size 
analysis shows that the gap {Eg) scales as Eg ~ h^^^^ 
where v{/^ = 0.25) = 1.18 ± 0.01^^^^^^ It is the scal- 
ing behavior of the second excited state which becomes 
gapped in the thermodynamic limit {N oo). Our data 
confirms the exponential decay of the specific heat for 
very low temperatures as < /ic- 

However, the general feature — which is the opening of 
energy gap due to breaking of rotational U{1) symmetry 
— has been appeared as the exponential decay of spe- 
cific heat at low temperatures. The exponential decay is 
vanishing as ^ /ic- At the critical point {h = he) the 
gap vanishes due to the presence of Goldestone modes 
which destroy the long-range antiferromagnetic order in 
^/-direction. For h > he the system enters the param- 
agnetic phase and the gap is proportional to the h — he 
which is clearly shown in Fig. O 

Moreover, the density plot in Fig. [2] shows three differ- 
ent regions at finite temperature which are distinguished 
by different energy scales. (I) The low temperature and 
h < he where the classical domain walls quasi particles 
define the dynamic of system^^. (II) The quantum criti- 
cal region where the dynamic is dictated by the quantum 
fiuctuations which become long ranged as both the tem- 
poral and spatial correlation length diverge. (Ill) The 
low temperature and h > he where the classical spin 
fiipped quasi particles represent the dynamical behavior 
of the model. 

We would like here to comment on the finite size effects 
in our calculations. Since the LTLM approach is based 
on several sampling in the Hilbert space, the finite size ef- 
fect will be weak and negligible for nonzero temperature. 
The size dependence of the specific heat is shown in Fig. 
© where different chains with TV = 10, 12, 14, 16, 18, 20 
have been considered. Unless the height of specific heat 
which is not the same for different in the low and high 
temperature regions all data fall on each other which jus- 
tifies the weak finite size effect. We ignore the small size 
dependence and claim that the our data of the specific 
heat can be interpreted for a very large chain length. The 
inset of Fig. (|3|) shows the small oscillations in very low 
temperature which is the result of level crossing where 
the ground state expectation is dominant in the parti- 
tion function and will vanish as the chain size becomes 
large. 




FIG. 2: The specific heat of XXZ chain in the transverse field 
(h) versus temperature (T) for A = 0.25. The data are from 
LTLM results and N — 20. (a) The three dimensional plot, 
(b) the density plot. (Color online) 



A remark is in order here. The finite size effect is weak 
and is compensated by random sampling if the model is 
not close to a quantum critical point. This is the case of 
our data presented in Fig. (|3|) at = 1. However, close 
to quantum critical points a finite size analysis should be 
implemented to find the correct behavior. For instance, 
the level crossing of the first two excited states should 
be incorporated in a finite size treatment as done in Ref. 
12 to find the scaling of energy gap. The situation is 
more complex to find the scaling behavior close to he 
since the value of he itself is obtained by a finite accuracy 
which is an extra source of error in determing the critical 
exponents. 

We have also calculated the specific heat of the XY 
chain, A = 0, which shows similar qualitative behavior 
to Fig. [21 It is in agreement with the quantum renormal- 



FIG. 3: The specific heat of XXZ chain versus temperature 
(T) in the transverse field h — 1 and A = for different chain 
lengths N = 10, 12, 14, 16, 18, 20. The inset shows the very 
low temperature regime. (Color online) 



FIG. 4: The density plot of ^/-component spin structure factor 
at momentum p — n versus the transverse magnetic field {h) 
and temperature (T). The chain length \s N — 20 and A = 
0.25. (Color online) 



ization group results^ ^ which states that the universality 
class of < A < 1 is the same as A = case. However, 
the critical field (he) is slightly smaller for A = which 
is approximately, hc{^ = 0) 3.1. 

C. Static Structure Factors 

We have calculated the static structure factor at finite 
temperature. The a-component spin structure factor at 
momentum p is defined by 

N 

G--{p) = ^{s?s?+,)e'^^ (9) 

X = l 

where (...) is the thermal average defined in Eq. (|6]). 

The Hamiltonian (Eq. ([1])) is gapless at = where 
there is no long range ordering in the ground state while 
the correlation functions decay algebraically, it is called 
spin fluid state. The transverse field {h 7^ 0) breaks the 
U(l) rotational symmetry to a lower Ising-like which de- 
velops a nonzero energy gap. The ground state then has 
long range antiferromagnetic order for < A < 1. How- 
ever, due to nonzero projection of the magnetization in 
the direction of the transverse field, it is a spin-flop phase. 

The antiferromagnetic ordering in ^-direction for the 
spin-flop phase implies that the ^/-component structure 
factor at p = TT diverges as the size of system goes to 
infinity. This can be seen as a peak in the structure 
factor. 

The density plot of G^^(7r) is shown in Fig. (j3j) versus 
temperature (T) and transverse field (h) for A = 0.25. 
The negative peak in G^^ (tt) (which is a signature of an- 
tiferromagnetic ordering) is the dark (blue) area in Fig. 
dH). It shows that the ordering is absent at = and 
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FIG. 5: The density plot of x-component spin structure factor 
at momentum p = versus the transverse magnetic field (h) 
and temperature (T). The chain length is A/" = 20 and A = 
0.25. (Color online) 



will gradually appears for increasing h. The antiferro- 
magnetic ordering in ^/-direction can be represented by 
staggered magnetization in the same direction {srriy) as 
the order parameter. Fig.(|4|) shows that srriy becomes 
maximum at h 2. The point h 2 can also represent 
the position where the gap becomes maximum. Further 
increasing the transverse field suppress the antiferromag- 
netic ordering until reaching the critical point h = he- 
It is clearly shown in Fig.(|4j) that srriy becomes zero at 
h = he. 

The density plot of x-component structure factor at 
p = has been presented in Fig.Q versus T and h. The 
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FIG. 6: The RPA three dimensional x^^(^) versus temper- 
ature for different transverse field, A = 0.25 and J± :^ 
0.00458 J for CS2C0CI4. For h < he the susceptibility di- 
verges at finite temperature while for h > he it only diverges 
at T = 0. (Color online) 



FIG. 7: The phase diagram of the bulk material in T (temper- 
ature) and h (magnetic field) plane. It is the result of LTLM 
for one dimensional chain and random phase approximation 
(RPA) to get the bulk property with inter-chain exchange cou- 
pling J± 0.00458 J for CS2C0CI4. The cross (x) represents 
the zero temperature critical field. 



value of G^^(O) is small for ^ 2 which implies low align- 
ment of the spins in the direction of the external field. 
However, this value starts to grow for h > 2 and obvi- 
ously saturates for h > he = 3.3. The maximum value 
of G^^(O) is the signature of long range order in the field 
direction where the moments saturates at the maximum 
value in the paramagnetic phase. However, as the tem- 
perature increases the ordering will be washed out due to 
thermal fluctuations which destroy the long range order 
in one dimension. The long range order at finite tem- 
perature is meaningful in the sense of three dimensional 
model which is composed of weakly coupled chains. This 
will be explained in the next sections via the random 
phase approximation (RPA). 

D. Magnetic Susceptibility 

Magnetic susceptibility can be calculated using, 

X?D = ^G""(fc). (10) 

To get the susceptibility of the bulk material which is 
composed of coupled chains with inter-chain exchange, 
Jx, we have used the random phase approximation 
(RPA). It is a mean field approach where the one di- 
mensional chain is treated exactly while the interactions 
between chains are weak and considered in a mean field 
approach. In the disordered phase the dynamical suscep- 
tibility of coupled chains at zero frequency is^^ 

^yy (j.) 

where Xid{p) is the one dimensional susceptibility in y- 
direction and momentum p which comes from the LTLM 



computations. We have plotted x^^iP = '^r) versus tem- 
perature in Fig. (|6]) for different transverse field. The 
susceptibility diverges at finite (nonzero) temperature for 
h < he which shows the phase transition from the disor- 
dered phase to the antiferromagnetic long range ordered 
phase. For h > hc^ the susceptibility only diverges at 
T = which justifies no long range ordered at finite tem- 
perature. 

We have also used the RPA to calculate the inter- 
chain exchange coupling {J±) from our numerical com- 
putations and experimental data^. This can help us 
to justify the mean field treatment of the interactions 
between chains. We have used the h = and Tn — 
2ijmK ^ 0.0813 J data presented in Ref. i@ and solved 
1 — J^Xid{Tn-,^) = to find J^. The result of this 
calculation is J±_ 0.00458J which can be compared 
with J±_ 0.0147J estimated in Ref. T^. The cou- 
pling between chains is two order of magnitude smaller 
than the intrachain interaction which verifies the RPA 
implemented here. We have then used the value of 
Jx ^ 0.00458J to find the whole phase diagram of the 
bulk material using our numerical data at finite magnetic 
field. The phase diagram of the three dimensional system 
has been presented in Fig. ([7j). The phase diagram shows 
the border between the ordered and disordered phase in 
the T — h plane. This is in good agreement with the 
experimental data presented in Fig. 11 of Ref. [5|. 

E. Summary and discussion 

We have implemented the symmetric algorithm of the 
Lanczos method to find the finite temperature (thermo- 
dynamic) properties of the anisotropic Heisenberg model 
in the presence of a transverse magnetic field. It is ar- 
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FIG. 8: The effect of number of random samples in X- 
component of spin structure factor at momentum p = tt versus 

temperate T-f^ /Tlnlnr nnlinp'i 
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FIG. 9: Specific heat vs. temperature with different number 
of random samples, with /i = and A = 0.25. The solid black 
line is the result of exact diagonalization. (Color online) 



gued that the symmetric algorithm^^ (LTLM) is more 
accurate at low temperatures than the finite tempera- 
ture Lanczos method^^. We have found that the specific 
heat decays exponentially at low temperatures for h < he 
which justifies the presence of finite energy gap in the 
lowest part of spectrum. The gap vanishes both at = 
and h = he with a maximum around h = 2. A similar 
behavior has been observed for the y-component struc- 
ture factor at the antiferromagnetic wave vector p = tt. 
It shows that the onset of transverse field opens a gap 
which stabilizes the antiferromagnetic (AF) order in y- 
direction. The AF ordering becomes maximum at h = 2 
and start to decrease by further increasing of the mag- 



netic field. The long range magnetic order vanishes at the 
critical point {h = he) where the gap becomes zero and 
the Goldestone mode destroy the ordering. The data of 
the one dimensional model can be used within a random 
phase approximation to find the phase diagram of a the 
experiments^ done on CS2C0CI4. We have estimated the 
interchain couphng to be J± o:^ 0.00458 J which justifies 
the RPA method. Moreover, we have plotted the phase 
diagram of the three dimensional model (weakly coupled 
chains) in Fig.© which is in a very good agreement with 
the experimental results presented in Ref,5. 

We have shown that the LTLM results are not finite 
size dependent if they are not close to the quantum crit- 
ical points. It is the result of random sampling which 
exists in the algorithm. However, a finite size analysis 
is required to find the scaling behavior close to quan- 
tum critical points. The scaling property of the energy 
gap has already been discussed in Ref.^, while there are 
other aspects which are open for further investigations. 

We would like also to comment on the number of ran- 
dom sampling in this method. We have plotted the x- 
component structure factor at p = tt versus temperature 
in FigE] for different sampling, R = 10, 20, 30 and the 
exact diagonalization on the full spectrum. It shows that 

= 30 is enough to get a good accuracy for the finite 
temperature behavior. A similar data for the specific 
heat versus temperature and R = 10, 20, 40, 60 in Figl9] 
shows that the number of R = 60 sampling reproduces 
well the exact results. We then conclude that our results 
is reliable for R = 100. One should note that the CPU 
required time is proportional to the number of sampling 
(R). 
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